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The current through ferromagnetic single-electron transis- 
tors (set's) is considered. Using path integrals the linear 
response conductance is formulated as a function of the tun- 
nel conductance vs. quantum conductance and the tempera- 
ture vs. Coulomb charging energy. The magneto-resistance of 
ferromagnet-normal metal-ferromagnet (F-N-F) SET's is al- 
most independent of the Coulomb charging energy and is only 
reduced when the transport dwell time is longer than the spin- 
flip relaxation time. In all-ferromagnetic (F-F-F) SET's with 
negligible spin-flip relaxation time the magneto-resistance is 
calculated analytically at high temperatures and numerically 
at low temperatures. The F-F-F magneto-resistance is en- 
hanced by higher order tunneling processes at low temper- 
atures in the 'off state when the induced charges vanishes. 
In contrast, in the 'on' state near resonance the magneto- 
resistance ratio is a non-monotonic function of the inverse 
temperature. 



I. INTRODUCTION 

The spin of the clcctjpn can provide new functional- 
ity in electronic devicesEJ and spin-polarized transport is 
therefore an active research helcL Spin-transport was pio- 
neered by Tedrow and Meserveyo and the discovery of the 
giant magnetoresistance effect (GMR) in metallic multi- 
layers has motivated a substantial rc-ncwed interest in 
this field. Recently spiuppolarized current injection in 
all-semiconductor devicesB and carbon nanotubesQ have 
been realized. The fabrication of tunnel junctions made 
of two ferromagnetic leads separated|-.lay an insulating 
layer has become well under controLufl These techno- 
logical advancements can make it possible to design new 
generation of electronic devices such as magnetic RAM's, 
sensors and ultimately perhaps quantum computers. 

In this article, we will focus on an elementary struc- 
ture, the ferromagnetic single-electron transistor and 
discuss the role of the ferromagnetic electrodes, spin- 
accumulation and Coulomb charging effects on the trans- 
port properties. A schematic picture of the system under 
consideration is shown in Fig. |l|. A normal (ferromag- 
netic) metal island is coupled to two ferromagnetic reser- 
voirs by tunnel junctions. There is a capacative coupling 
C between the metal island and the metal reservoirs. The 



island is so small so that the charging energy Ec = e^/2C 
associated with the addition of a single electron to the 
island can be larger than both the temperature and the 
applied source-drain bias. The electrostatic potential on 
the island can be controlled by a gate voltage Vg ca- 
pacitively coupled via a capacitance Cg to the island. 
We consider an island weakly coupled to the reservoirs 
by many channels with transmission probabilities much 
less than unity, but the total tunnel conductance G can 
be smaller or larger than the quantum conductance Gk- 
In the simplest case, the single-electron transistor (SET) 
comprising of tunnel junctions with junction resistances 
much larger than the quantum resistance Rk — h/e^ ex- 
hibits Coulomb blockade of the single electron tunneling 
at low temperatures when the gate voltage vanishes, and 
a conductance peak when CgK) = e/2. For spin-polarized 
charge transfer the different transport properties of the 
spin-up and spin-down electrons have to be taken into ac- 
count. Therefore the behavior of the SET is determined 
by the interplay between i^jpin-dependent transport and 
Coulomb charging effects. Q^Ej 

All- ferromagnetic (F-F-FJl SET's have recentijubeen 
investigated experimentallyQ and theoretically.l3^ElEj The 
tunneling magneto-resistance (TMR) in F-F-F systems 
is defined as the relative conductance difference by 
changing the magnetization of the ferromagnetic island 
from being parallel to anti-parallel to the magnetiza- 
tions of the ferromagnetic reservoirs (which are kept 
parallel). Ferromagnetic islands have a short spin- flip 
relaxation time and several works have computed the 
transport prouer|ties of F-F-F systems disregarding spin- 
accumulation.QElllj An enhanced TMR in the Coulomb 
blockade regime was found experimentally,B and it was 
sub-sequentially shown that the TMR could be doubled 
due to co-tunneling@ and that higher order tunneling prfc 
cesses could give even larger enhancements of the TMR.E^I 
In what follows, we will use a non-perturbative approach 
to calculate the tunneling magneto-resistance (TMR) of 
the F-F-F SET showing that even larger enhancements of 
the TMR than what has been reported in Ref . |l3| is possi- 
ble. We will also discuss the TMR in the 'on'-state when 
the Coulomb blockade effect is lifted. It will be demon- 
strated that the TMR in this regime is a non-monotonic 
function of the temperature. 

The physics of the TMR in F-N-F systems is com- 
pletely different from the mechanism of the TMR in F- 
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F-F systems. F-F-F systems exhibit a finite TMR even 
in the absence of any spin-accumulation on the ferromag- 
netic island. The transport properties can be modeled by 
an all-normal metal device with magnetization configu- 
ration dependent tunnel conductances. In contrast, the 
TMR of F-N-F systems is uniquely related to the spin- 
accumulation on the island. Spin-accumulation governs 
the transport properties and leads to a measurable TMR 
when the spin-flipj-iielaxation time Tgf is longer than the 
dwell time so thatli3 



nf/{xh)>GK/G, 



(1) 



where x is the spin-susceptibility, G is the tunnel con- 
ductance and Gk = e^/h is the quantum conduc- 
tance. In order to realize spin-accumulation a small spin- 
susceptibility is required (^) which implies a small nor- 
mal metal island with a large Coulomb charging energy 
Ec and consequentULjCoulomb blockade effects should be 
taken into account .IliJ When the tunnel conductance G is 
much smaller than the quantum conductance Gk and 
the temperature is not too low, transport can be de- 
scribed by sequential tunneling processes and this was 
performed in Refs. p[-pl]. However, Eq. (|l|) shows that 
spin-accumulation is easier to realize when the tunnel 
conductance is large, G > Gk- In this strong tunneling 
limit the quantum fluctuations of the charge on the is- 
land as well as the spin-acccumulation should be taken 
into account with a theoretical description that goes be- 
yond the low order perturbative treatment presented in 
Refs. pl p^ The tunneling magneto-resistance (TMR) in 
F-N-F systems is defined in a different way than for (F- 
F-F) systems: The TMR is the relative conductance dif- 
ference on going from a configuration where the magne- 
tizations in the ferromagnetic reservoirs are parallel to a 
configuration where the magnetizations are anti-parallel 
(ffp ~ 5ap)/<7p- Changing the magnetization configura- 
tion in F-N-F systems only changes the spin-dependence 
of the tunnel conductances but conserves the total tun- 
nel conductances. Hence, a much weaker coupling to the 
charge degrees of freedom of the island and a reduced in- 
fluence of Coulomb charging should be expected. In the 
Coulomb blockade regime at low temperatures the total 
conductance through the system is suppressed. Conse- 
quently, also the transport of spins into the island de- 
creases and at sufficiently low temperatures the dwell 
time is longer than the spin-flip relaxation time so that 
(|l|) is no longer satisfied. The spin-accumulation and 
thus the TMR must therefore be reduced in the Coulomb 
blockade regime which is contrary to the behavior of the 
TMR in F-F-F systems. We will show that the linear 
response TMR in F-N-F systems in the Coulomb block- 
ade regime equals or is smaller than the classical high- 
temperature TMR, which is in contrast to the results for 
F-F-F systems. 

We present in this paper a unified analytical expres- 
sion for the linear response conductance of F-F-F and 
F-N-F set's with dirty metallic islands. Our formula is 



valid to all orders of the ratio between the tunnel con- 
ductances and the quantum conductance, for arbitrary 
spin-flip relaxation times and takes into account general 
band-structures for the ferromagnetic reservoirs and the 
normal/ferromagnetic metal island. Our results are in 
this sense general. The analytical expression for the lin- 
ear response conductance is used to discuss the magneto- 
resistance of F-F-F and F-N-F SET's. 

The paper is organized in the following way. The sys- 
tem is described and the model Hamiltonian introduced 
in Section |l| In Section III the linear response conduc- 



tance is calculated, including the multi-band effect of the 
electrodes and the spin-accumulation on the central is- 
land. The TMR is evaluated in Section [V in the case 



of ferromagnet-normal metal-ferromagnet SET's and in 
Section ^ in the case of all-ferromagnet SET's. Section 
VI concludes the paper. 



II. MODEL 

We consider a ferromagnetic single-electron transistor 
comprising of a normal/ferromagnetic metal island con- 
nected to ferromagnetic leads by tunnel junctions. The 
central island is capacitively coupled to a gate voltage Vg 
via a capacitance Cg. The Hamiltonian of such a device 
is 



H= (Hi + H, + Hr]+Hc+( Hi, + H„ + h.c 



(2) 



The quasi-particles in the left reservoir are described by 
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(3) 



where n and m are the band-indices, cr is the electron 
spin, 1 is the momentum in the left reservoir, W{^™ deter- 
mine the quasi-particle energy bands and the hybridiza- 
tion between the bands for a given spin, and cj'^ creates 
an electron with spin cr and momentum 1 in band n in 
the left reservoir. The hat Q denotes an operator. There 
are similar Hamiltonians for the quasi-particles on the is- 
land {I — > i) and the quasi-particles in the right reservoir 
(Z — > r). The Coulomb charging effects are included via 



H, 



(4) 



where the induced electron number is riox = GgVg/e and 
the spin-dependent number of excess electrons on the is- 
land is NX = yii„ c^^cV-^. We will below also make use of 
the spin-dependent excess number of electrons in the left 



(right) reservoir, = ^m^il^^a iK = Er„ ^ 



rn ra ra 



The tunneling between the reservoirs and the island 
is taken into account via the tunneling Hamiltonians 
Hpi = J2a ^pic {P = ^ denotes left and p = r denotes 
right) : 
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(5) 



The tunneling matrix elements allow tunneling be- 
tween different bands with the same spin. Spin-flip scat- 
terings during the tunneling processes are disregarded, 
but can be included in the formalism leading to a renor- 
malization of the spin-asymmetries of the junction resis- 
tances. 

We consider the situation when the temperature is 
much larger than the level spacing and the level spac- 
ing is much smaller than the Coulomb charging energy. 
The level spacing is consequently taken to be continous. 
Co-tunneling and higher order tunneling ppj-ocesses have 
both elastic and inelastic contributions ,EZl e.g. in the 
case of co-tunncling the elastic contribution means that 
the same electron tunnels through both of the two junc- 
tions whereas the inelastic contribution corresponds to 
two different electrons that tunnel in the two junctions: 
One tunnels into the island above its Fermi level, and 
another jumps out of the island from below the Fermi 
level. The elastic contribution is sensitive to the phase 
of the electron on the island and consequently becomes 
vanishingly small for dirty systems when the character- 
istic time of tunneling through the macroscopic barrier 
caused by the Coulomb energy Ti/Ec {Ec = e^/(2C)) is 
larger than the classical time L'^/D for diffusion through 
the island (L is the. size of the island and D is the dif- 
fusion coefficient ).E3 We consider the inelastic transport 
regime exclusively and disregard the elastic contributions 
to the tunneling processes, i.e. for an estimate of a dif- 
fusion constant on the island oi D — lOcm^/s and a 
charging energy of -Ec = 1-?^ we assume that the size of 
the island L < O.lfim (in reality there is also scattering 
at the boundary of the island that further relax this as- 
sumption). The occupation of the energy levels on the 
island can then be described by an energy e and spin a 
dependent non-equilibrium distribution /^(e). The non- 
equilibrium distribution /^(e) is spin-dependent to allow 
a non-equilibrium spin-accumulation and only equals the 
equilibrium Fermi-Dirac distribution f{e — fi) when the 
bias voltage vanishes {fi is the equlibrium chemical poten- 
tial). Furthermore, the leads and the island are metallic 
so that the tunnel conductances (see Eq. ( |l6| ) below) are 
energy independent on the scale of the Coulomb charg- 
ing energy which is much smaller than the Fermi energy. 
We consider the linear response regime where the bias 
voltage is the smallest energy scale in the system. Under 
these assumptions, only the total occupation of spin-up 
and spin-down electrons on the island, J de[f^ — f{e — fi)], 
determine the macroscopic Coulomb barrier and the to- 
tal effective tunneling rates. Consequently to all orders 
in perturbation theory of the tunneling Hamiltonian, we 
can describe the electron occupation on the island by a lo- 
cal chemical potential /i^ = /i+ / de[/^(e) — /(e — /i)] since 
the energy-dependence of the distribution f^{e) does not 
influence the transport properties. It is thus e.g. ir- 
relevant if there is substantial energy-relaxation on the 



island so that the non-equilibrium distribution /^(e) ap- 
proaches the Fermi-Dirac distribution /(e— /i^) or if there 
is no energy-relaxation so that the non-equilibrium dis- 
tribution /^(e) is a linear combination of the Fermi-Dirac 
distributions in the left and right reservoirs with different 
local chemical potentials. 

The electrons in the right reservoir and in the left 
reservoir are in local equilibrium with chemical poten- 
tials ^\t) and /z''(t), respectively and as seen above the 
electrons on the island can be described as in local equi- 
librium with a spin-dependent chemical potential /^^(i). 
The temperature is the same in all subsystems. The ther- 
mal average of the spin-dependent current from the left 
{p — I) or right {p — r) reservoir to the island can then 
be written as 

im = zeTr [(4,.(t) - i?^,(t)) eM~(3K{t))\ (6) 

with the non-equilibrium time-dependent grand canoni- 
cal potential 



K{t) = H- y^iii'mi + f,''^{t)K + f^iit)K 



(7) 



The spin-dependent chemical potentials on the island /i^ 
and /ij^ should be determined from the flux of particles 
aiid spins into the island. The spin-flux into the island 



(/j + /[)-(/! + /[) = 



Tsf 



Current conservation through the system requires 
{I\+I[) + {I]- + II)=0. 



(8) 



(9) 



The excess number of spins on the island s is related to 
the non-equilibrium difference in the chemical ppteiitials 
by the spin-susceptibihty Xs via s — Xsi^J■\ ~ fJ-'l)!^^ For 
non-interacting electrons Xs = D, where D is the density 
of states. The expression for the current (m and the 
conservation laws for spins (^ and particles (^j uniquely 
determine the spin-dependent local chemical potentials 
on the island, /ii^ and fii^, and consequently the current 
through the device. We will focus on the linear response 
regime in what follows. 



III. LINEAR RESPONSE CONDUCTANCE 

In the linear response regime a perturbation expansion 
in terms of the small differences in the chemical poten- 
tials on the left reservoir, on the island and on the right 
reservoir can be performed: /i'(t) — IJ- + 5^^{t), fJ-^{t) = 
fi + Sf/" (t) , and /i^ (t) = fj, + Sfj-l^ (t) . The non-equilibrium 
grand canonical potential is K{t) = Kq + SK{t), where 
the time-independent unperturbed grand canonical po- 
tential is = -ff - fJ-T^A^a + + K) and the 
time-dependent perturbation is SK = —J2ai^f^''^cT + 
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Sfi'^N^ + Sfil^N^). We perform a time-dependent uni- 
tary transformation with the unitary matrix U{t) = 
exp[i J_^dti5K{ti)]. The transformed grand canonical 

potential becomes Kt{t) = Kq + 6H{t), where to the 
lowest order in the non-equilibrium chemical potential 
differences 



SH{t) ^-Y^lKt) f dh [Sfi'^ih) - Sf^lih)] 



(10) 



and the sum over the parameter q corresponds to the left 
reservoir q — I and the right reservoir q = r. The spin- 
currents through the left and the right barriers can be 
found from a standard linear response calculation: 



mu;) = A^n^«(c.) [6^,^u;)-S^,lic.)] , (11) 
euj ^-^ 



qs 



where the retarded current-current correlation function 
is 



(12) 



This correlation function will be evaluated using path 
integrals. The computation is similar to the case of ul- 
trasmall tunnel junptipBS or metallic SET that has been 
studied extensivelyJla^tS The main steps in _the_ calcula- 
tion of the correlation function is as follows.li3£3 We use 
the equilibrium finite-temperature Matsubara formalism 
to find the retarded current-current correlation function. 
The imaginary time generating functional is defined as 



Z[{V}] 



D[{c*,c}]D[^]e- 



(13) 



with the generalized action 



Sg = - dTTTpncrC*{T)drC{T) 

Jo 

+ / drKoir)- f dTY,Il[T)ijl[T). (14) 

A Hubbard-Stratonovich transformation is performed to 
reduce the quadratic term in the charging energy to a 
linear coupling-of the particle number operator with an 
auxiliary field.E£l After integrating out the electronic de- 
grees of freedom and including the principal contribution 
of the tunneling processes the generalized effective action 
is derived. The current-current correlation function can 
then be found by performing the functional derivative 
nP|(T,T') = Z-^d"^ Z / {driP{T)dri'i{T'))\r,^o. The result- 
ing Fourier transform of the current-current correlation 
function is 



-nil - ^ [5,„.qs9l90 + . 

The dimensionless tunnel conductance (p = /, r) is 



(15) 



n' 



''pa 



(16) 



Here Df^ , D^^ and Df^ are the spin and band dependent 
densities of states and the renormalized transmission co- 
efficient at the Fermi energy is ipa- = (C/po-)^ tpa-Uia, where 
Upa (Uia-) is the material-specific matrix that diagonal- 
izes the band Hamiltonian in reservoir p (island) for elec- 
trons with spin a. In a single-band model the conduc- 
tance is simply proportional to the density of states in 
the reservoir and the island, gP — Dp„Di„\ta\'^ ■ The re- 
lation is in general more complicated as shown in (p^). 
The correlation function gQ^iUn) is 

go = — / dre-"-x(r)r(r) (17) 
«w„ J 

with the phase-phase correlation function 



1 

A;— — oo 

D[ip]e-^^'^^ cos [(^(t) - ip{0)] . (18) 

bk 

The effective action is 

l2 



S[ip] = / dT 



Ml 



g-' / dr / dr'x (r - r') cos [^(r) - ^(r')] . (19) 



with the boundary condition — ¥^(0) + 27rfc. 

Since the leads and island are (ferro-)metallic with con- 
tinuous spectra, the damping kernel xi^) is an even 
function of the imaginary time t with a period f3 and 
the Matsubara- Fourier components xi^i) = — |wi|/47r. 
The coupling of the island to the reservoirs is deter- 
mined by the algebraic total dimensionless classical con- 
ductance g'^' = X^po-ffS- The partition function is Z = 
J2klb -^[v] 6^P(~'S'M)- Another correlation function in 
(H) is 



9. = ^f dre-^^^ X 

ILUn . Z 



E 



2?Me-5[^l/t(^,T)/t(</,,0) (20) 



where 



It{ip, ^) = y drx (ct - t) sin Yp{a) 



(21) 



Due to the current conservation (^ the correlation func- 
tion gi (|2^) does not appear in the final form of the 
stationary current. 

Knowing the current-current correlation function, we 
can find the non-equilibrium chemical potential on the 
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island by using the conservation laws for spins 
particles (^. We find the spin-accumulation 

/^T ~ _ {9\9l- g\9\)/{9]gi) 
and the average chemical potential 



.9 M 



and 



(22) 



(23) 



.9T + 5i ^ 

where the spin-flip conductance (jsf is related to the spifc 
susceptibility Xs and the spin-flip relaxation time Tgt viaE2l 
ffsf = T^Xs/Tsf- From the non-equilibrium chemical poten- 
tials on the island we obtain the current through the sys- 
tems by using the relation for the current (O) and finally 
we calculate the conductance g = el/{5ir— Sf/'): The 
linear conductance in units of the quantum conductance 
e^/h of the FSET is 



9 = 



9^9' 
9'' +9' 



790 



I - 



(glfff ~ 9\9l) /{g^'g 9]9i) 
1 + 5sf [(5o5t)"^ + {gogi)~ 



(24) 



where and {g^ and are the spin-up and spin- 
down conductances of the left (right) junction in units 
of the quantum conductance e^/h, and we have also 
introduced the combined conductances g' = g' + g', 

5*^ = 5t + gl' g^ ^ g\ + g^ and .gj. = 5x + 9l- The re- 
tarded correlation function go ( p^ ) is the same as for an 
all-normal metal single-electron transistor. It renormal- 
izes the conductance due to the Coulomb charging effects 
and depends on the induced electron number CgVg/e, 
the ratio of the charging energy to the temperature PEc, 
and the coupling between the island and the reservoirs 
through the total conductance g^ + gK Through Eq. ( |2^ ) 
we have thus reduced the problem of finding the con- 
ductance of ferromagnetic single-electron transistors to a 
well-studied problem of the renormalized conductance in 
all-normal metallic single-electron transistors. The sec- 
ond term in the bracket of (|4|) is a signature of the spin- 
accumulation and vanishes when Tgf — > {gst —>■ 00). 
The spin-accumulation causes a reduction of the con- 
ductance. The resulting conductance (^) can be un- 
derstood in terms of the equivalent circuit shown in Fig. 
(^, which is identicai-to the high-temperature classical 
circuit for two spins jl3 but now generalized to include 
Coulomb charging effects by renormalizing the conduct- 
naces so that g\ gog\, g\ -> 50.9', 9\ -> 909] and 
9\ ~* go9\- The spin- flip conductance is not renor- 
malized, since spin-flip processes do not change the num- 
ber of electrons on the island and consequently not the 
charging energy. We will now use (|2j) to find the TMR 
ratio in F-N-F and F-F-F SET's. 



IV. F-N-F MAGNETO-RESISTANCE 

We first focus on F-N-F systems. The quantity of in- 
terest is the TMR [g-p — gAp)/gp- When the magnetiza- 



tions in the left and right leads are parallel, the tunnel 
conductances can be characterized as gl = g^l + sP/)/2 
and gl = g''{l + sP^)/2, where s = +1 (s = -1) de- 
notes spin-up (spin-down) and Pi — (g' — g\)/{g\ +5') 
and Pr = [g^ — g\)/{g\ + g\) denote the polarizations of 
the left and right tunnel junctions, respectively. In the 
anti-parallel configuration the tunnel conductances are 
g[ = c,'(H-sP/)/2andg^ = g^'il- sPr)/2. The total con- 
ductances of each junction, = g\+g\ and g^ — g\ + g\j 
are thus independent of the magnetization configuration. 
Hence, the correlation function go and consequently also 
the prefactor in (|2^ ) remain invariant on going from the 
parallel to the anti-parallel configuration. The magneto- 
resistive effect is thus solely due to the spin-accumulation 
which is contained in the second term of Eq. (p4|). The 
TMR of F-N-F systems is 



gp - gAP 
gp 



pp. 



(25) 



where the parameter = 4(7sf, 
reduces the spin-accumulation and hence the TMR due 
to spin-flip relaxation, 7^ = (g' — g^)/{g^ + g^) (0 < 
7g < 1) reduces the TMR due to the asymmetries in the 
conductances of the left and right junctions, and (3g = 

{P^l + Ig? + P.^(l - Ig? - 2PiPr{l - ll))lA (0 < Pi < 
1). The TMR is proportional to the polarization of the 
left and the right tunnel conductance. This result (|25| ) 
reduced to the high temperature or large bias result in 
Ref . |l^ when go ^ 1. 

The spin-flip relaxation time is governed by spin-orbit 
coupling and magnetic impurities at low temperatures. 
We assume that Tgf is temperature-independent below 
the Coulomb charging energy (at very low temperatures 
this assumption can become invalid, but that is a regime 
beyond the investigation in this article). Firstly, since 
the correlation function go can only decrease due to 
Coulomb charging effects the factor ag appearing in ( |25| ) 
can only increase at low temperatures and the TMR in 
the Coulomb blockade regime is equal to or smaller than 
the classical high-temperature TMR. Secondly, when the 
spin-flip relaxation time is much larger than the dwell 
time, the TMR is independent of the charging energy, 
temperature and tunnel resistances for practically all 
temperatures since ag — > except at extremely low tem- 
peratures. Thirdly, in this regime we expect the nonlin- 
ear response TMR to have a re-entrant behavior: The 
linear response TMR is identical to the TMR when the 
bias is much larger than the Coulomb charging energy 
which equals the classical high-temperature TMR. In the 
intermediate tuas-cegime effects of the Coulomb charging 
energy appearQ'EHlj in the sequential tunneling limit and 
an enhancement of the TMR has been demonstrated in 
Ref. using the co-tunneling formalism. Finally ( |2^ ) 
and (E5[) can be generalized to finite frequencies by let- 
ting go go{^) and 1/Tsf — *■ l/rgt + ito so that the AC 
response can be used to detect the spin-accumulation.E2l 
Our general result (pa) is consistent with the calculations 
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performed in the sequential tunneling regime 




V. F-F-F MAGNETO-RESISTANCE 

The spin-flip relaxation time in a ferromagnet is_tjyp- 
ically much shorter than the transport dwell timeJlJ In 
the limiL^Tst — > the conductance of F-F-F systems re- 
duces toEj 

I r 

g=4^9oif3E,,g' + g\n,^). (26) 

The TMR of F-F-F systems is defined as 

^ ^ ^AP - _ gp - ffAP ^27) 
-Rp 5AP 

where Rp = RK/gp {Rap = Rx/gPiP) is the resistance 
of the FSET and gp (gxp) is the dimensionless conduc- 
tance of the FSET when the magnetizations in the leads 
and the central island are parallel (anti-parallel). (Note 
that we use a slightly different definition for the relative 
tunneling magnetoresistance in F-N-F systems than in 
F-F-F systems in order to make our discussions coherent 
with previous work). For F-F-F systems it is assumed 
that the magnetizations in the leads remain parallel, and 
that only the magnetization direction of the central is- 
land changes by applying an external magnetic field. 

The TMR is caused by the magnetization configu- 
ration dependence of the total junction conductances 
g^ — g\ + g[ and g^ — g^ + g^^ and therefore also by 
the magnetization configuration dependence of the renor- 
malization factor go . The latter dependence gives rise to 
the enhancement of the TMR in the Coulomb blockade 
regime. The first term of the effective action (^9|) govern- 
ing the renormalization of the conductance go (O) only 
depends on the Coulomb charging energy of the device, 
and thus in the weak tunneling limit go — > 1 and the 
dimensionless conductances and the TMR of the FSET 
reduce to the classical values. The second term in the ac- 
tion (p^ ) is proportional to the sum of the dimensionless 
conductances of the tunnel junctions and consequently is 
magnetization configuration dependent. 

In order to evaluate the TMR the conductances of the 
FSET in different alignments must be calculated. We 
characterize the polarizations of the junctions by the rel- 
ative polarizations P and P' associated with the reser- 
voirs and the island, respectively. In a simple two-band 
model with spin-independent tunneling probabilities the 
polarizations P and P' are directly related to the density 
of states in the reservoirs and the island. In general there 
is not such a one-to-one correspondence between the po- 
larizations of the tunnel barriers and the polarizations 
of the density of states, but we can still characterize the 
spin-dependent asymmetries in the tunnel conductances 
by the parameters P and P'. We denote the dimension- 
less conductance of the left (right) junction averaged over 
the parallel and anti-parallel alignments by g^^^ . In the 



parallel configuration the spin-dependent conductances 
are 

gZ = l9'^''i^ + ^m + sP'), (28) 

where s — 1 (s = —1) for spin-up (spin-down) electrons. 
Hence, the total classical dimensionless junction conduc- 
tance is 

9p=gp + 9p = g{i + PP'), (29) 

where g = <f + g^ ■ Similarly, in the anti-parallel config- 
uration the conductance is 

g%.s = \g'''{i + sP){i-sP'). (30) 

The total classical dimensionless junction conductance is 
thus 

gtp^9ip+9lp^9{^-PP')- (31) 

In the high temperature limit, one can substitute the 
forms of gp and g^Jp described above into the definition 
of the TMR (p7|), and obtain a simple expression 7^^' = 
2PP'/{1 - PF)iov the classical TMR. 

If the thermal energy is larger than the Coulomb charg- 
ing energy, the path integrals can be evaluated semiclas- 
sically, and the lowest order quantum corrections to the 
TMR can be cJaljOiJied. In this regime the conductance 
of the FSET isElN 

5o(w = 0) ~ 1 - ^ + (0.0667 0.0185<7,f)(/?£:c)^ 

-^e--^/'^^-</2cos(2^ne.) (32) 

with 77 = P{AP) for the parallel (anti-parallel) align- 
ment. Substituting this expression into the definition of 
the TMR (|^), we obtain 

7- = 7^1 + (1 + 7^1) [pE^f + , (33) 

where = 0.02(7p, and the leading term of the gate 
voltage dependent part of the semiclassical TMR is 

57(ne.) = |^(l+7^')e-'^^/'^^^x 

[e-9Ap/2 _ g-sp /2] cos(27rnex) ■ (34) 

In Ref. ^ riox = was considered and consequently 
the gate voltage dependent contribution (^ ) was disre- 
garded. For large junction conductances or high temper- 
atures, this term is exponentially suppressed and the in- 
fluence of the gate voltage on the TMR is negligible. The 
gate voltage dependence part of the TMR is important at 
low temperatures. Although Eq. ( |3^ ) is not accurate at 
very low temperatures, it can be used to estimate when 
<57(7T-ex) becomes relevant. For FSET's with positive po- 
larizations PP' > the conductance in the anti-parallel 
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configuration is smaller than the conductance in the par- 
allel configuration g'^p < gp, and the maximum TMR 
occurs when ricx = 0. The TMR attains its minimum at 
n-cx = 1/2, and returns to the maximum at riox = 1 since 
the TMR is a periodic function of with period 1. 

At low temperatures, an analytical expression of the 
TMR is not available. We will study the maximum and 
minimum value of the TMR at low temperatures via 
Monte-Carlo simulations. We use Ni leads with P = 0.23, 
and Co island with P' = 0.35 in all numerical calcula- 
tions below. 



A. TMR in the 'ofF'-state 

In the absence of induced charges, the geometric phase 
factor 2Trikncx in the phase-phase correlation function 
( p^ vanishes. The density matrix for each winding num- 
ber k is then positive and the phase-phase correlation 
function can be computed directly using standapl Monte 
Carlo simulation techniques. In our earlier worktj a large 
enhancement of the TMR, roughly 4 times larger than 
the classical TMR for g = 12 and Ec/ksT = 40, was 
found. In this section we will demonstrate that the TMR 
is further enhanced at even lower temperatures. We will 
also address whether even larger junction conductances 
increase or decrease the TMR ratio. 

The enhancement of the TMR ratio can be clearly 
demonstrated in the co-tunneling regime, .wbp.re the con- 
ductance can be calculated analytically.EZlEfl Since co- 
tunneling is a second order tunneling process and the 
first order sequential tunneling vanishes in the Coulomb 
blockade regime, the conductance of the FSET is pro- 
portional to the square of the classical conductance with 
a temperature dependent prefactor. Consequently the 
TMR (|2^) is enhanced |-to a value slightly larger than 
twice the classical TMRB A naive extension of this result 
is that beyond the sequential and co-tunneling regimes, 
the TMR is enhanced by a factor of n when the n-th 
order tunneling process dominates. However this picture 
is too simplified since contributions from different orders 
of tunneling have different temperature-dependent pref- 
actors and one cannot in general specify the tunneling 
process that dominates for given junction resistances and 
temperature. In the strong tunneling regime perturba- 
tion theory breaks down and a number of higher order 
tunneling processes have to be taken into account simul- 
taneously. 

When the thermal energy is not much smaller than 
the charging energy, the semiclassical formula (^) can be 
modified to fit the TMR values by replacing the Coulomb 
charging energy by a renormalized charging energy. This 
substitution agrees well with the reaiilts of earlier investi- 
gations on the single-electron box.EElEj The renormaliza- 
tion scheme also predicts that the effective charging en- 
ergy at sufficiently low temperature decreases for larger 
junction conductances due to the coupling to the leads. 



This means that for a given temperature, the TMR in 
the strong tunneling regime is smaller for larger junction 
condiictances. The role of strong tunneling is therefore 
(ifia/.EJ The constructive role is to increase the TMR by 
the inclusion of higher order tunneling processes at low 
temperature. The destructive role is caused by the renor- 
malized charging energy to a smaller value E* <C Ec so 
that the TMR is closer to the classical value at a given 
temperature. 

When the thermal energy is much smaller than the 
Coulomb charging energy strong tunneling is manifested 
as an enhancement of the TMR by lowering the tem- 
perature. We compute the TMR values for two values 
of the classical total dimensionless conductance g — 12 
and g = 24. For a given conductance g and temperature 
the phase-phase correlation function (nSl) is computed via 
Monte-Carlo simulations for a sufficiently large number 
of discrete points in the imaginary time interval [0, /?]. 
The phase-phase correlation functions after two succeed- 
ing N samplings are compared. If the curves are smooth 
and sufficiently close to each other, they predict roughly 
the same conductance after the analytical continuation 
of the Matsubara- Fourier components. The small differ- 
ence determines the uncertainty of the final result, which 
is shown as the size of the data point in the figures below. 
When the deviation between the two phase-phase correla- 
tion functions is larger than the requested accuracy they 
are averaged to give the phase-phase correlation function 
for the case of 2A^ samplings, and a further simulation is 
performed for the next case of 2A^ samplings. This pro- 
cedure is followed until the final result converges. Using 
the computed phase-phase correlation functions, the con- 
ductances of the TMR can be calculated via the Kubo's 
formula Eq. (|^. The analytical continuation from the 
Matsubara frequency to th&*eal frequencyE3 is performed 
via the Fade approximate.^ 

The resulting TMR ratio in the absence of the gate 
voltage is shown in Fig. (||) for g = 12 and 24. In both 
cases, the TMR values at low temperature are larger 
than the classical values, which shows that the TMR 
is enhanced by higher order tunneling processes in the 
Coulomb blockade regime. For ^ = 12 at ksT/Ec — 0.01, 
the TMR is enhanced by a factor of 8 compared to the 
classical value, which is almost twice as large as the value 
that we reported for ksT/Ec = 0.025 in our previous 
work.Ej For a given temperature, the TMR for g = 12 
is larger than the one for 5 = 24 in the temperature re- 
gion shown in Fig. (^), which means that the smearing 
of the Coulomb charging energy via higher order tunnel- 
ing processes is stronger for larger junction conductances. 
Therefore for experimentally easily accessible tempera- 
tures, only moderate enhancement of the TMR for the 
FSET formed by tunnel junctions with extremely large 
conductances can be obtained. In order to experimen- 
tally observe large TMR ratios in SET's tunnel conduc- 
tances of the order of 10 times the quantum conductance 
are recommended. 
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B. TMR in the 'on'-state 

In the 'on'-state when ricx — CgVg/e — 1/2 module 1, 
there is no Coulomb blockade effect even when the ther- 
mal energy is much smaller than the charging energy. 
In the sequential tunneling regime, the conductance of 
the FSET approaches halt-of the classical value when 
the temperature is lowered.cJ Beyond the sequential tun- 
neling regime, the reduction of the conductance of the 
FSET by lowering the temperature is slow at relatively 
high temperatures, but at sufficiently low temperature 
the conductance attains values smaller the conduc- 
tance in the sequential tunneling regime .BEj The conduc- 
tance in the 'on'-state in the strong tunneling limit and 
its implication for the TMR will be investigated in this 
section. 

In order to calculate the 'on'-state conductance in the 
strong tunneling regime by Monte-Carlo simulations, we 
reformulate the phase-phase correlation function as 

r°"(r) =r'=(T)/Z", (35) 

where the function r''°(T) is the mean value of the phase 
factor 

r°(T) = — V / Dipe-^'^^^ cos(7rfc) x 

cos[^{t) - vp(O)] , (36) 

with respect to the partition function which is a di- 
rect summation of the density matrices for all winding 
numbers fc, Z° — YlkL^oo Ibk Dfe~^^'^^ and the parti- 
tion function at resonance is renormalized in a similar 
way, 

^" = ^ E / Dipe-s^^^ cosink) . (37) 

In the above formulas, both the numerator and the de- 
nominator of the phase-phase correlation function r°'^(T) 
can be calculated with respect to the density matrices of 
the FSET with a geometric factor which equals one. The 
computation can therefore be performed in the same way 
as in the case of vanishing gate voltage. The complica- 
tion is that we now have to calculate both r''°(r) and 
instead of a single correlation function. 

When the thermal energy is not much smaller than 
the charging energy, quantum fluctuations are small and 
the simulations are rather fast. Consequently the stable 
phase-phase correlation functions can be readily obtained 
using the approach described above. A typical result is 
shown in Fig. (jj). At lower temperatures, the imaginary 
time P is longer, and more lattice points are needed to 
guarantee the convergence of the Trotter product. More- 
over, the computed correlation functions for the same 
amount of sampling numbers turn out to fluctuate more 
strongly at lower temperatures as shown in Fig. (a). For 



5 = 12 at f3Ec > 60 and g = 24 at jSE^ > 100, the 
correlation functions computed via the approach of dou- 
bling the sampling number still fluctuate considerably 
within a realistic computing time. Therefore we do not 
attempt to perform simulations for temperatures lower 
than (3Ec — 100. At the lowest temperature the correla- 
tion function at the points with large fluctuations is re- 
placed by a best fit smooth function so that the analytical 
continuations in the Kubo's formula lead to a converging 
conductance of the FSET. The accuracy of this smooth- 
ing approach needs to be determined by long time Monte 
Carlo simulations in the future, since we cannot guaran- 
tee that the single data point at the lowest temperature 
is accurate within its size. 

The TMR ratio as a function of the inverse dimension- 
less temperature is shown in Fig. (^. For g — 12, the 
TMR first increases by lowering the temperature. The 
values are almost the same as the case of the vanishing 
gate voltage. Around PEc — 10, the gate voltage depen- 
dent contributions to the TMR becomes important shift- 
ing the TMR value for CgVg = e/2 farther away from 
the value for CgVg = at lower temperatures. There- 
fore while the TMR in the Coulomb blockade regime is 
enhanced further at lower temperatures, the TMR at res- 
onance decreases rapidly to the values smaller than the 
classical TMR for f3Ec > 30. Negative TMR values could 
possibly occur at even lower temperatures, which needs 
to be determined quantitatively by further investigations 
which at present cannot be performed due to the com- 
puting time required. For g = 24, in contrast to the 
case of smaller junction conductances g — 12, the TMR 
increases constantly and slowly by lowering the temper- 
ature down to a value of O.OlEc/ks- The curve more 
or less follows the tendency of the TMR in the case of 
vanishing gate voltage. This is understandable because 
the gate voltage dependent contributions to the TMR is 
proportional to the factor exp{—gfl/2) as expressed in 
(p^), which is smaller for larger g, therefore the differ- 
ence between the maximal TMR at CgVg = and the 
minimal TMR at CgVg = e/2 determined by the gate 
voltage dependent contribution is also smaller, and the 
variation of the minimal TMR will change the direction 
at lower temperatures. In both cases, the difference be- 
tween the maximal TMR and the minimal TMR is found 
to be larger at lower temperatures, therefore a wider vari- 
ation range of the TMR can be achieved by tuning the 
gate voltage. 

VI. CONCLUSION 

We have derived a general formula for the linear con- 
ductance of single-electron transistors containing ferro- 
magnetic elements. In F-N-F structures, the TMR is 
almost independent of the Coulomb charging energy and 
is only reduced at sufficiently low temperatures when 
the effective transport dwell time becomes shorter than 
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the spin-flip relaxation time of the normal-metallic is- 
land. In F-F-F systems, we have calculated the magneto- 
resistance as a function of gate voltages for a wide range 
of temperatures in the strong tunneling regime. In the 
'off'-state, the magneto-resistance can be enhanced at 
low temperatures by higher order tunneling processes. 
However, higher order tunneling processes also reduce the 
effective charging energy, and slow down the rate of the 
enhancement of the magneto-resistance when the junc- 
tion conductances are much larger than the quantum con- 
ductance. Consequently the largest magneto-resistance 
for a given temperature is obtained when the junction 
conductances are larger, but not very much larger than 
the quantum conductance. The magneto-resistance of 
the device in the 'on' state has a more complicated tem- 
perature dependence. Due to the contributions from the 
induced charges, the magneto-resistance increases slowly 
for very large junction conductances. In contrast to this, 
the magneto-resistance is reduced well below the classi- 
cal magneto-resistance for not very large junction con- 
ductances at sufficiently low temperatures. 
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FIG. 1. The single-electron transistor. The central island 
is connected to the left and the right reservoirs by tunnel 
junctions. The electrostatic energy of the island can be tuned 
by a gate voltage capacitively coupled to the island. 



FIG. 2. The equivalent circuit of the low-bias conductance 
of the FSET including the spin accumulation effect of the 
central island. 



FIG. 3. Tunneling magneto-resistance for CgVg = as a 
function of the inverse of the dimensionless temperatures in 
the strong tunneling regime for g = 12 (filled dotted), and 
for g = 24 (filled diamonds). The errors are smaller than the 
sizes of the symbols. 
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FIG. 4. Phase-phase correlation function as a function of 
the imaginary time for g = 12, l3Ec = 20 and CgVg = e/2. 
The dotted Uncs are the Monte-Carlo data for 10^ samplings, 
each. The dot-dashed line is the average of the dotted lines, 
the dashed line is the Monte-Carlo data for next 2 x 10® sam- 
plings, and the solid line is the average of the dot-dashed line 
and the dashed line. 

FIG. 5. Phase-phase correlation function as a function of 
the imaginary time for g = 12, pEc = 60 and CgVg = e/2. 
The dotted lines are the Monte-Carlo data for 10® samplings, 
each. The dot-dashod line is the average of the dotted lines, 
the dashed line is the Monte-Carlo data for the next 2 x 10® 
samplings, and the solid line is the average of the dot-dashed 
line and the dashed line. The smooth long dashed line in 
the figure is used to replace the fluctuating solid line to cal- 
culate the Matsubara-Fourier components of the correlation 
functions. 

FIG. 6. Tunneling magneto- resistance for CgVg = e/2 as a 
function of the inverse of the dimensionloss temperatures in 
the strong tunneling regime for g — 12 (filled dotted), and 
for g = 24: (filled diamonds). The errors are smaller than the 
sizes of the symbols except for the point at lowest temperature 
for each set of parameters, where the error estimation is not 
accurate. 
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